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Abstract. It is shown that a lagrangian system whose Legendre 
transformation degenerates along a hypersurface behaves in a strange 
manner by jumping from time to time without any "visible cause" . 
In such a jump the system changes instantaneously its coordinates 
as well as its momenta. The mathematical dscription of the phe- 
nomenon is based on the theory of impact, refraction and reflection 
developed by one of the authors and the observation that a hamil- 
tonian vector field, understood as a relative one, can be associated 
with any lagrangian, degenerated or not. Necessary elements of 
the general theory of such systems are reported and a detailed de- 
scription of a post-relativistic oscillator showing such a behaviour 
is given. 
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1. Introduction. 

In Lagrangian mechanics two, in a sense, extremal situations were 
widely studied. One of them, classical, corresponds to systems with 
nowhere degenerated Legendre map. In this situation the Legendre 
map identifies the Lagrangian dynamics with the corresponding Hamil- 
tonian one. On the contrary, in the second case the Legendre trans- 
formation is supposed to be everywhere degenerated (and of constant 
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rank). Such a situation is common in gauge theories and by this rea- 
son was studied in several works, starting with the pioneering paper by 
Dirac ]I[. But in a generic, in the sense of singularity theory, situation 
the Legendre map L is almost everywhere non-degenerated except a 
critical hypersurface S. Since C maps the tangent bundle T (M) of 
the configuration space M of the system in question to the cotangent 
bundle T* (M) , which both have the same dimension, the most generic 
singularities of C are just folding points (see, for instance, ||). In such 
a case S subdivides T (M) or a suitable domain of it into a number of 
regular regions, inside each of which the dynamics admits the standard 
description either in Lagrangian or in Hamiltonian form. But what 
happens when a trajectory reaches the singular hypersurface SI In 
this paper we show that it jumps. More exactly, if its arrival point is 
of the fold type , i.e. a generic one, then the moving point in M at 
a certain instant disappears suddenly and at the same time reappears 
in one or more, generally, distant points to continue its smooth mo- 
tion up to an eventual subsequent jump. It is worth to emphasize that 
collisions and impacts of any kind in mechanics as well as refractions 
and reflections in geometrical optics in the context of the hamiltonian 
approach are phenomena of the same mathematical nature. In this 
case the "jumping" quantity is the velocity, but not the space-time 
coordinates (see 0, appendix 3). 

In this paper, which is the first of a series planned to be written 
in this connection, we present the general mathematical background 
and then illustrate it with an example of a certain physical flavour, a 
relativistic oscillator in the post-Galilean approximation. Our approach 
is a fusion of two simple but, it seems, important observations.. The 
first of them is that the hamiltonian vector field can be associated 
naturally with any lagrangian, independently of the fact whether the 
corresponding Legendre map is degenerated or not. This field, however, 
is a relative one (see below) with respect to the Legendre map. The 
second observation is an analogy with the hamiltonian theory of impact 
as developed by one of the authors ( ||, appendix 3). 
It was not our goal here to discuss variational aspects of the problem. 
It will be done in a separate purely mathematical paper. We only 
remark that variational problems with non- regular lagrangians were 
already studied by K. Weierstrass and various authors later on (see, 
for instance, the Weierstrass- Erdmann theorem in 0). 

It seems that not too much attention was paid to physical systems the 
actions of which are not regular everywhere. We mention here recent 
works by G. Vilasi, I. Pavlotsky and their collaborators (j7|, [|J), where 
some particular results concerning the post-Galilean oscillator and the 
Darwin two electrons-model were obtained. Our approach is, however, 
completely different and leads directly to a complete dynamical picture. 
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We hope that the forthcoming study of more realistic models will 
clarify the physical content of the proposed mechanism. At the moment 
its eventual applications look so attractive that it would be wiser to 
postpone speculations. 



2. Hamiltonian Theory of Impact: the Transition 

Principle. 

Here we describe with the necessary details discontinuous hamilton- 
ian systems, which will help us to understand the behaviour of dynam- 
ical systems described by singular lagrangians. Namely, the geometry 
that describes jumps of phase trajectories for such hamiltonian systems 
is the same as for lagrangian ones with singularities of fold type. In 
other words, the principle controlling these "jumps" in the pure hamil- 
tonian case is quite transparent and shows what one has to do in the 
more complicated lagrangian situation. 

Let ($, Q) be the phase space of a dynamical system with Q = 
^2i dpi A dqi being a symplectic 2- form on $. Suppose then that $ 
is divided by a hypersurface T into two closed domains $+,$_, having 
T as their common boundary, i.e. <9$ + = T = <9$_. Suppose also 
that the hamiltonian of the system is smooth on $±. In other words, 
if H-t = H\$> ± , then H± G C7°° ($±). This, in particular, means that 
H±\t G C°° (r), but it is not supposed that H + \r coincides with if_|r- 
So H, and consequently the hamiltonian field Xh associated to it, is 
well defined on $\r and is bi-valued on I\ 

As it is well known, the local coordinate expression of X H± is 



^dH± d dH ± d 



dpi dqi dqi dpi 

The corresponding canonical equations : 

• _ dH± _ _dH± 

Qi — K j Pi Fj ) 

opi dqi 

describe the motion of the system inside $±.But when the phase trajec- 
tory arrives at V it must " decide" under control of which hamiltonian to 
proceed on. The transition principle ( || j prescribes how this decision 
should be taken. Below we recall this "recipe". 

First, remember that to each point x G V the characteristic line 
lx C T x ($) of T at x is associated ([§). This one- dimensional subspace 
is the skew-orthogonal complement of the hyperplane T X (T) C T x (<&) 
with respect to the bilinear, skew-symmetric, non-degenerate form Q x . 
In other words: 
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Figure 1. In- and out-points 



Since l x C T X (T), for any iGT, this way we get a one- dimensional dis- 
tribution x i— > l x on T , whose integral curves are called characteristics 

of r. 

Example. Let S be a hypersurface in a manifold M, dimM = n. 
Then , denoting by it : T* (M) — > M the canonical projection, we 
obtain the hypersurface T = n^ 1 (5) C $ = T*(M). Its characteristics 
are straight lines contained in the fibers T* (M), q G S. In fact, let 
q n = be the equation of S in a certain local chart. Then this is also 
the equation of T in the corresponding chart (q,p) on $. By definition 
a vector £ G T x . ($) , x 6 T, is collinear to the characteristic direction 
at x iff covector fi x (£, •)vanishes on T x . (r), i.e iff ^ (£, •) = \dq n for 
some A G R. It is easy to see that vector -S-\ x satisfies this condition. 

Hence characteristics are the integral curves of the vector field ^-, i.e. 
straight lines. 



Let x G T. We say that x is a +-in-point (resp., a +- out-point) 
if Xh + \x is directed toward $ + (resp. $_). Similarly, we say that 
x is a — in-point (resp., a — out-point) if X^_| x is directed toward 
$_(resp. $ + ) (see Fig. |l|). In- and out-points of if± are separated by 
a hypersurface r^cT along which X H± is tangent to T. 

Suppose now that the phase trajectory, starting from a point inside 
of $ + (resp. $_), reaches a point x G T at an instant t, and let E 
be the constant value of H + (resp. iJ_) along the phase trajectory for 
t <t . Denote by j x the characteristic curve of T passing through x, 
and by S^ (resp. E^) the hypersurface {H + = E} (resp., {iJ_ = E}) 
of $ + (resp., $_). A point y G j x r\ Sj (resp. 7 x .flEg) is called decisive 
for a; if it is an +- in-point (resp. a — in-point). Now we can state the 
following 
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Figure 2. The Transition Principle (y and z are deci- 
sive for x) 



Transition Principle: When a moving point reaches the separation 
hypersurface T at a point x with energy E, it continues then its mo- 
tion from all x-decisive points y G E £ simultaneously under control of 
corresponding Hamiltonians H±. The passage of the phase point from 
x to the y's is assumed to be instantaneous. 

The transition principle is illustrated in Fig. 0. 

In || some examples illustrating this principle in mechanics and in 
geometrical optics are given. To these we will now add the following 
one. 

Example (laws of reflection and refraction). 

Let M be an inhomogeneous, anisotropic optical medium, referred 
to a system of orthogonal coordinates (qx,q2, 93)- Denote by V(q) = 
c/n (q) the light velocity at point q. Here c is the light velocity in 
the vacuum and n(q) stands for the refraction index. Then (|J) the 
propagation of the light rays is described by the canonical system: 



with 



(1) 



qi— H Pi , Pi— —H q 



H(q,p) = V(q) 



i = l,2,3 
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where ||p|| = JJ2 Pi- 



Let now M + ,M_ be two isotropic optical media with refraction in- 
dexes n+(q) and n_(g) respectively, separated by a surface S. Sup- 
pose that a ray starting from a point in M + reaches S at a point 
q = (qi,q2,%) with an impulse p = (Pi,p 2 )P3)- The corresponding 
velocity is u = (t7i, ^Us), with ^ =^= H Pi (q,p) = V {q)Pi/ 



Pose: M = M + U M_ , $ = T* (M) , $± = T* (M±) and T = U 

T* (M). So, the hypersurface T separates $ + from $_ and the phase 
trajectory of the ray reaches T at a point x = {q,p). By the previous 
example the characteristic curve 7^ of S passing through x is a straight 
line contained in the fiber T± (M). Let us choose coordinates (gi, g 2 , 93) 
in such a way that S is tangent at q to the hyperplane {g 3 = q 3 } and 
the g3-axis is directed toward M_. Then the characteristic direction at 
x is parallel to -£- and 7^ is described by equations 

Qi(t) = Qi 
(2) ft(t) = pi + 6 a t , ieR, z = l,2,3 

If E = H + (q,p), then, in view of (|2|) and (d), the intersection 7^flS^ 
is composed of two points x and x* corresponding to values t — and 
t* = — 2p 3 , respectively The point x* = (q,p*) with p* = (pi,p 2 , — P3) 
is decisive for x. In fact, by the choice of coordinates: 

(X H+ ) x , (q 3 ) = ^± (x*) = -V + {q)^ = ~ (**+)*(*) < ° 



Therefore, the transition principle tells that the reflected ray does al- 
ways exist, and starts from the same point q G S in the direction 
v* = V (q) p* / \\p* \\ corresponding to p*. Further, v, v* are coplanar 
with the normal to S at q and form with it equal angles (j), ip + Respec- 
tively In fact 

\ v t\ \Pl\ IPs I / 

cos0 = -i-^-r = -f-^- = ^-4r = cost/m 

\\v \\ \\P \\ 



On the other hand the intersection 7^nS £ corresponds to the values 
of t satisfying equation 

(3) t 2 + 2p 3 t + (l-n 2 ) \\f\\ =0 

with n = n-(q)/n+(q). Equation ([3]) admits real solutions iff ( & 
f-) > 1 - n 2 , i.e. iff 



^3 
sin (j) — Wl : — 2 < n 



^+ (?)' 
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Therefore, for values of <fi greater than arcsinn there is no refraction 
("total reflection"). If instead < arcsinn (or if n > 1), then (|^) has 

two real solutions, t = —p 3 ± p 3 , with p 3 = \\p\\ * / ( & + n 2 — 1 j , to 

which correspond the two intersections (q, Pi,p 2 iPa), (q,P~x,P~2i ~ Pa) -Of 
these, only the first one is decisive for x (one can see it, as before, by 
considering the orientation of Xjj_). Therefore, for < <fi < arcsinn 
the refracted ray does exist, and its (initial) direction p = (jPi,p~2,pa) 
is coplanar with the incident ray and the normal to S at q, and forms 
with this an angle ip- suc h that (Snellius' law): 



A. 

sm< _ 

n 



sin ip_ _ l\- M. 



3. Lagrangians with Singular Hypersurfaces. 

In this section the Hamiltonian theory of impact and refraction, 
as described in the previous section, is extended to the singular la- 
grangians. This will allow us, as it was already mentioned in the intro- 
duction, to describe discontinuities of the motion that occur when the 
phase point of the system reaches the singular surface. We start with 
recalling the concept of relative vector field. 

3.1. Relative vector fields. Let M and iV be differentiable man- 
ifolds, related one another by smooth mapping F : M — > N. An 
R-linear operator X : C°°(iV") — ► C°°(M) is called a relative vector 
field on iV along F if it satisfies the Leibniz rule: 

X(fg) = F*(f)X(g) + F*(g)X(f) f, g G C°°(N) 

If / £ C°° (M) and X is a relative vector field, then fX is also a such 
one. Therefore relative vector fields on iV along F form a C°° (M)- 
module. Denote it by T>(N,M;F). An "absolute" vector field X on 
A^ can be considered as a relative one along the identity map id^ . 

Geometrically, a relative vector field X ET> (A, M; F) can be thought 
as a map which associates to any point a G M a vector X a tangent to 
A at the point F (a). Namely 

X a (f) ■= [X(f)\ (a) , / G L7°°(A) 

Conversely, to each differentiable map a : M — > T(N) such that 
p N o a = F ,where p^ '■ T(N) — > A is the canonical projection, one can 
associate the relative vector field X defined by: 

[X(f)](a):=a(a)(f) , a G M, / G L7°°(A) 

If x = (xi, , x m ), y = (y±, , y n ) are local coordinates on M and 

A respectively, the corresponding local expression of X G T>(N, M; F) 
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is: 



n ft 



i=l 



where functions X l (x) are the components of X x with respect to 

the basis -^- \f{x)i ,^~\f(x) of T F ^(N). Below we will use the 

simplified notation 

having in mind that the vector field (|]) acts on a function <j) (y) G 
C°° (N) as follows: 

X(0)(xi,....,x m ) = J2x i (x) — (y 1 (x), ,y n (x)) , 

with yi = yi (x) being the coordinate expression of F. 
For more details concerning relative vector fields, see ||. 

3.2. Singular Lagrangians and Relative Hamiltonian Theory of 
Impact and Refraction. Let M, dimM = n, be the configuration 
space of a dynamical system described by a lagrangian L G C°°(T(M)), 
and let £ : T(M) — > r*(M) be the corresponding Legendre mapping. 
Recall that, in a fixed local chart (gi, ...., q n ) on M, C is represented by 
equations: 

(5) qi = qi , i = l,....n 

Pi = L Vt (q,v) , i = l,....n 

where (q,v),(q,p) are special coordinates on T(M) and T*(M), respec- 
tively, associated with (gi, ,<?n)- 

Let S 1 be the singular points locus of the Legendre map : 

S :={xe T(M) | rk d x C < n} 

This hypersurface can be described as follows. The Jacobian matrix of 
C at point (q, v) is: 



(6) d{q,v)£ 



1 



where L vq = \\L Viq \\ and L vv = \\L ViV \\ is the hessian matrix . Since 
the determinant of (|B|) is TC := det L vv , the equation of S is: 

W(g,v) = 

We will assume that S is a regular hypersurface of T(M). This assump- 
tion implies that for any point (q, v) G S, d( q , v )Tt does not vanish. 
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According to the standard procedure, the motion of the system out- 
side S can be described by Euler-Lagrange equations: 

(7) J * = Vi 

i — 1, , n. Equations (^) can be rewritten in the normal form: 

(8) Vi=fi(q,v) 

in a neighbourhood of any point (q, v) G T(M)\S. On the other 
hand, this is no longer possible if (q, v) 6 S. Namely, accelerations are 
undetermined on S and velocities may have discontinuities. 

Recall (see, for instance, 0) that x &T (M) is a fold point of C if 

(9) Ker d x C n T x (5) = {0} \/x e S 

(10) 4W ^ OVieS 

Below we restrict our analysis only to fold points, which for a generic 
L form an open everywhere dense subset in S. 

In the next section we will see that the lagrangian of the post-galilean 
oscillator exhibits also some irregular singularities, which can be stud- 
ied as well. 

For our purposes it is important to note that the hamiltonian vec- 
tor field associated with a lagrangian exists even if the corresponding 
Legendre map possesses some singularities. But in such a case this 
hamiltonian field becomes a relative one. Namely, define the relative 
vector field X r H el e D (T* (M) ,T(M); £) by: 



- W -|- + £M,»)|- 



;ii) *»' = £».£: + £Ms, 



If £ is (locally) regular, i.e. a (local) diffeomorphism, then 

x H = (c- i yox% 1 

is the standard hamiltonian field associated with L.This is easily seen 
from the standard coordinate expression of a hamiltonian field: 

In fact, in this case 

(12) H(q,p) = E(q,v(q,p)) , 

where v = v(q,p) is the inverse of the Legendre map and 

E(q,v) = ^2viL Vi (q,v)-L(q,v) 



V 



J 



10 F. PUGLIESE, A. M. VINOGRADOV 

is the energy. Keeping in mind @ one finds 

Hence, the vector field Zl corresponding to Euler-Lagrange equations 
for L is well-defined on T (M) \S and, locally, 



^=z(^|-+^)|-) 



with /j's given by (§). If [/ C T (M) \S is such that £[/ d = C\ v : £/ -»• 
£ (LQ is a diffeomorphism, then the image X v G £> (£([/)) of Zx,|j/ 
with respect to Ly is well defined: 

Xu = {C v l ) o Z L o C* n 

and, moreover, is a hamiltonian vector field (in C(U)) corresponding 
to the Hamilton function H v = E o C~ l e C°° (£ (U)). If J7' is another 
regular domain for C such that C {U') = C (U) , then, generally, Hyi ^ 
Hu or, equivalently, Xu ^ X V i (in C{U)). In other words, in C(U) 
a multi-valued hamiltonian field is well-defined. So, on the whole, a 
multi-valued hamiltonian field is defined and its various branches are 
matching one another along £ (S). 

Now we pass to describe the analogues of in- and out-points (see 
section 2) in the considered context. The following elementary facts 
from singularity theory (see, for instance, (@)) are needed for this 
purpose. First, recall the notion of a submersion with folds. Let P, Q 
be two manifolds, with dimP = n > dimQ = m. The first- order jet 
[F] of a map F : P — > Q at a point x G P may be considered as a triple 
(x, y,p), with y — F (x) and p : T x (P) — ► T y (Q) being the differential 
of F at x. The manifold of all the first-order jets of maps from P to Q is 
denoted by J 1 (P, Q). The submanifold Si C J 1 (P, Q) is composed of 
all triples (x, y,p) such that rkp = n — 1. The map j^P : P —>■ J 1 (P, Q) 
sends a point x G P to [P] x . P is called a submersion with folds if: 

1) j 1 F is transversal to Si; 

2) .fTer d^P is transversal to Si(F) = (j 1 F)~ 1 (Si), for every x G 

It is easy to see that in the case n = m condition 1) is equivalent to 
the fact that the jacobian of P has only simple zeroes along Si(F). 

The local structure of a submersion with folds is described by the 
following 

Theorem 3.1. Let F : P — > Q be a submersion with folds and let 
xG Si(F). Then there exist coordinates (xi, ....,x n ) on P and (yi,....,y m ) 
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on Q, centered at x and F(x) respectively, in terms of which F takes 
the form: 

?/i = a?i 



Vm—l %m—l 

— " I 2 I I 2 

Dm % m ^ x m+l ^ ±x n 



< see [[J ►. 
In the case we are interested in, P = T(M), Q = T*(M), F = C, dim P = 
dimQ = 2n. Assumptions @, (pT0|) guarantee £ to be a submersion 
with folds. According to theorem j3.1| , its normal form is: 

V\ = xi 

(13) - 

2/2n-l — %2n-l 

_ 2 
V2n — X 2n 

It results easily from the above description and from ( [L3|) that the range 
of C locally belongs to the half-space yi n > . This allows us to extend 
the definition of in- and out-points and the corresponding transition 
principle to the lagrangian case. Namely, a point x G S is called an in- 
point if Xjj l \ x is directed toward the range of C,Xjf\ x (jfen) > , while 
it is called an out-point if X^f l \ x is directed outside of it, X^f l \ x (jfen) < 
0. 

Let us note that the pullback C* (Q) of the canonical symplectic form 
on T* (M) along the Legendre map is degenerated on S and is of rank 

def 

2n — 2 at any fold point. The restriction Q$ — £* (^) \s continues to be 
of rank In — 2 due to (0) . This means that the kernel l x of Qs at a fold 
point x G S is one-dimensional. This way one gets a one- dimensional 
distribution on Sf id- Characteristic curves are integral curves of it. 
Denote by 7 X the characteristic curve passing through x G S. 
Now it is clear how to extend the transition principle to the lagrangian 
case. Namely, calling decisive for x G S any in-point y G 7 X belonging 
to the same level Eg of energy of x, the principle can be stated as 
follows. 

Transition Principle (lagrangian case). When a phase point 
moving along Zl reaches at an instant a point x G S, it then continues 
its motion along all trajectories of Zl issuing from points decisive for 
x. Moreover, the passage from x to a decisive point is instantaneous. 

Remark 1. Note that the transition principle implies that the en- 
ergy does not change under an impact with S. 

Remark 2. According to the principle there are in general as many 
possible phase trajectories after the impact with S at a point x as are 
the points decisive for x. 
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4. Relativistic Oscillator 

In this section the behaviour of a relativistic oscillator whose la- 
grangian possesses fold singularities is analysed on the base of Transi- 
tion Principle. We have chosen this example to show the above theory 
in action mainly because of its relative simplicity: its phase trajectories 
and characteristics can be described analytically without difficulties. 
However, we will see that even in this simple case the behaviour of 
phase trajectories with respect to the singular surface is rather inter- 
esting, at least from a geometrical point of view. 

By applying the transition principle to describe the discontinuities 
of motion of relativistic oscillator we discover a rather remarkable phe- 
nomenon. Namely, if the energy exceeds a certain level and at the 
same time the velocity is not too high, the oscillator starts jumping. 
In other words, after a smooth motion it instantaneously changes its 
position (as well as velocity). In the classical example of reflection and 
refraction of light, or elastic collision of bodies and particles, there is 
no discontinuity of position. It is worth noting that the "jumping" 
motions of the oscillator are in a good consistency with the smooth 
ones. 

Finally we note that the relativistic oscillator possesses also singu- 
larities of non-fold type. These are very degenerated and the phase 
portrait in this region is rather curious. 

4.1. Relativistic Oscillator. Recall that there were proposed various 
relativistic generalizations of the standard harmonic oscillator (see for 
instance 0). Some of them possess singularities, while others not. 
Below we study the two-dimensional post-galilean oscillator of tensor 
rank 2 (J7|), possessing both fold and not fold type singularities: 

(14) L = L(r,x) = —mc 2 



/r^+l-J (i + 2 



x 



Here m, r are the mass and the characteristic length of the oscillator 
respectively, linked by the relation r$ = \/]jr c (with k being the elastic 
constant); r is the distance between the oscillating mass and the elastic 

2 

force centre; x = ^ is the square of oscillator velocity, measured with 
respect to the light velocity c. If we fix in the plane of motion a system 
of orthogonal coordinates ((71,92) with the origin at the centre of the 
force, then obviously: 



9 9 V \+ V 2 

q{ + qi x= — 

Note that 



(15) 0<x<l , 

due to the fact that v 2 < c 2 . Below we refer to M = R 2 = {(qi, ^2)} 
as the configuration space. So the lagrangian ([TJ|) is defined in the 
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domain U C T (M) = {(q, v)} = R 2 x R 2 , denned as: 

U = {(q,v) e T(M) \qeM, \\v\\ < c} 

However we will often not distinguish between 14 and T(M). A similar 
convention will be adopted also for the cotangent bundle T*(M)). 

Below we will systematically use in U the system of coordinates 
(r, 0, x, 0) (or equivalently (r, 0, x, u), with u = 9 — <fi), where <fi and $ 
are the angle between gi-axis and r = (51,^2), and the angle between 
r and the velocity vector v = (fi,t> 2 ), respectively. 

To simplify general considerations concerning lagrangian (|14|) it is 
convenient to work with a generic lagrangian of the form 

L = L (r, x) 



(16) 

The energy function 

(17) E = viL vl + v 2 L V2 
takes the following form for lagrangian 

(18) E (r, x) = 2xL x - L 
which in the case of oscillator fll~4] ) becomes: 



L 



E(r,x) 



mc 



VT 



+ 



X 



Lagrangian (^) admits also another integral of motion, namely the 
angular momentum: 



(19) I(r, x, u) = -zL x {r, x) {q x v 2 - q2Vi] 



I(r,x, u) 

For lagrangian flTJ]) it is specified as: 

1 



-\fxL x [ 



r, x sin u 



I(r, x, u) = mcryx 



VT 



X 



sinu 



The integral I corresponds, via Noether's theorem, to the infinitesimal 
symmetry : 

d_ 

dq 2 d(p 

of lagrangian (|16|) . 



X 



d 8 

oqi 



4.2. Singular hypersurface of relativistic oscillator. The Legendre 
map £ associated with lagrangian (0) is given by: 



(20) 



Qi ~- 


= Qi 




Pi ~- 


- L Vi - 


c z 



1,2 
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Figure 3. Singular surface (r = r/rO) 



The corresponding jacobian matrix (|6]) in terms of standard coordinates 
(q, v) and (q,p) in T (M) and T* (M) , respectively, has the entries: 



7 ^<?j 



c 2 r 



-Vtqj 



2 4 



i,j = l,2 



So, the corresponding Hessian is 



H{q,v) = L VlVl L 

It is easy to see that: 

(21) 

and, therefore, 

Hence, 



- T 2 

«2«2 ±J v 1 v 2 4 



-jL x (L x + 2xL xx ) = H (r, x) 



L x + 2xL xx — hi x 



H (r, x) = —L X E X 



E X L X 







is the equation of the singular hypersurface S. In other words 

s = s x u s 2 

with Si = {E x = 0},S2 = {L x = 0}. Each of these hypersurfaces is 
fibered in tori (0, 6), some of which may reduce to circles or to a point, 
depending on L. The bases of these "fibrations" are curves Ci,C 2 
in (x, r)-plane, given by equations E x (r, x) = and L x (r, x) = 0, 
respectively. These curves for the oscillator ([TJ) are shown in Fig |3|. 
In this case the intersection 7 = S\ fl 52 is the circle {r = r ,x = 0} 
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included in the null section M C T(M). 7 is the locus Sing S of 
singular points of S. This follows easily from the relations 

riq. = — ri r = — (L xr E x + L x E xr ) 

Vi Vi 

v i ~2 2 \ xx x ' J-'x-t-'xx) ) 

and from the fact that Si\y and 5*2 \7 are regular. 

Now we pass to describe how the kernel of the Legendre map behaves 
along the singular hypersurface S. It follows from (||) that a vector 

£ = X/»=i ( a *^" + ^^ ) belongs to the kernel of gLC iff 

(22) ai = a 2 = , L vv I * J = f Q 

In full details, the second of conditions f|22]) looks as: 



'1 

( 2 2 \ 2 

(23) I L z + -glraa-Vi I 61 + -jLaa-Vi^^ 

L xx v 1 v 2 b 1 + I L x + —L xx v\ ] 6 2 - 



o -—axe "± ^z w i 1 1 -"a; ' o 



By construction these equations are linearly dependent at any point of 
S. 

For a point (q, v) G S 2 \S\ there are two possibilities: L xx (q, v) = 
or L^ (q, v) 7^ 0. In the first case the system becomes trivial and 
Ker d( q , v )£ is two-dimensional. This never happens for the oscillator 
(|i4|). In the second case system (|23| ) reduces to the linear equation: 

vih + V 2 b 2 = 

Except for the points (v\ = v 2 = 0), Kerd( q>v )£ is one-dimensional. 

It is generated by the vector v 2 -J^ u i££T> which coincides with — ^ 

in coordinates (r,<p,x,u). This shows that in both cases Kerd( q>v )C 
is tangent to S 2 . So fold type singularities do not belong to S 2 and 
dim C (S 2 ) < 3. For the oscillator ( ^4|) C (S 2 ) is 2-dimensional and 
coincides with the null section of T* (M). 

Now we go to describe fold points belonging to Si\S 2 . These form 
an open domain S[° C Si everywhere dense in Si. 

If (q, v) G Si\S 2 , it follows from (|2T| ) that L x = —2xL xx , hence (p3|)2 
becomes 

v 2 bi - v\b 2 = 

Therefore we have that Ker d(q )V )£ = Span {v\-J^+v 2 -J^-) = Span ( J^) 
on Si\S 2 . 

Thus, singularities of C along S\ are of a substantially different na- 
ture from those along S 2 . Namely, the kernel of dC is transversal to 
Si\S 2 (due to the fact that E xx 7^ on it) and is tangent to S 2 , since 
^- = on it. Hence, dim C(S 2 ) < 3 and the transition principle 
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can not be applied to S2. In the next section we shall see that even 
characteristic directions are undetermined on S2. 

4.3. Characteristic curves on S. To simplify computations we will 
make use of coordinates (r, x, (f>, 1?) (or, equivalently, (r,x,u, </>)). Let 

P = ^Pidqi 

i 

be the universal 1-form on T*(M) (||). Then 

2 

(24) £*(p) = L Vl dqi + L V2 dq 2 = -\fxL x (r, x) (cos udr + r sin ud<p) 

c 

It follows from ([?§) that £*(p) | 5a = 0. Hence 

r(n)| Sa =o 

Therefore, £ (S2) is a lagrangian submanifold in T* (M) (with possible 
singularities) and characteristic directions on 5*2 are undetermined. 

Now we pass to describe characteristics on S(° . Since gradTi 7^ 
on it, the equation E x (r, x) = of Si can be solved with respect to one 
of the variables, say r: 

r = 7"i (a;) on Sj° 

In the case of the oscillator ( |I~4] ) : 

r i ( x ) — 3" 

(l-x)f 

So (x, m, 0) can be taken as local coordinates on S(° . Then from (f24l) 
we get: 

2 r- r ' 

£*{p) \sx= - v%L x (ri(x),x) r 1 (x) cos udx -\- rAx) sin u d> 
c L 

so that 

(25) C*(n)\ s = dC*(p)\ Sl = 

2 

- [a(x, u)dx A du + /3(rr, w)cfc A c?0 + 7(2;, -u)<iu A d<p] 



with 



\fxL x (r x \ 



x),x)r x \x) sin it 



(5(x,u) = — \r\{x)\fxL x (ri(x),x)l sinw 
dx 

7(2, w) = ri(x)y / i;L a ,(r 1 (a;), x) cosu 

Characteristic directions on S± are described by a characteristic vec- 
tor field X G £> (Si), i.e. such that 

(26) r(fi)| Sl (A»=0 
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If 

r\ r\ r\ 

ox ou o<p 



then (EH) is equivalent, in view of (pop, to 





This system is of rank 2 on S 1 / , and its fundamental solution is 

(a,b,c) = (7,-/5, a) 

Therefore, 

_ d d d 

dx du d(j) 

ft (1 (i 

q(x) cos it— q (x) sin-u— — \- r-, (x)\/xL x (ri(x) , x) sinw— — , 

OX OU 0(j) 

with q(x) = ri(x)yfxL x (ri(x),x) and characteristic curves are solu- 
tions of the system: 

x = q(x)cosu 

(27) it, = —q (2) sin u 

■ r[ (x) . 

<p = q(x) — — sinM 

7"i (X) 

Integration of system ( ^7|) is reduced, obviously, to its subsystem 
(p7|)i, 2 , whose solutions are to be described in the rectangle [0, l[x [— n, it] 
due to the ciclicity of u. 

It follows from (pj])i,2 that 

du ( x) 

(28) — = -^-ftanu 

dx q (x) 

and, consequently, 

/du , 1 • 1 f 1 i x ) , f dq , / \ 1 
= In sinw = — / — — — dx = — / — = — In \q (x) \ + const. 
tanw J q(x) J q 

= In rT ^r , c > 
\q(x)\ 

Hence the general integral of (p8|) is 

a 

(29) sinw = , , a E R 

g(x) 

For the oscillator ( [TS ) it is specified as: 

2a (1-z)* 



smti 



mc 2 r x | 
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Figure 4. Characteristic curves on SI 



nil 




Figure 5. Intersection between a chracteristic curve 
and the energy level surface 



Let us remark that q(x)sinu 



- I 

2 



So (1291) shows that I 



is 



constant along characteristic curves of Si. Therefore, by the transition 
principle, angular momentum (as well as energy) does not change after 
the impact with S\. 

The curves (^9]), denote them by j a , for the oscillator are shown 
in Fig. |j. Since the variable u is cyclic mod In and 7 a and 7_ a are 
symmetric with respect to the x-axis we can limit ourselves to dealing 
with the curves in the rectangle (x,u) G [0, l[x[0, 7r]. 

Let P = (x, u, 4>) G Si and let 7a be the characteristic passing 
through P, E = E (P) . The intersection between the energy level 
surface Ej and Si is the torus Tg = {r = n (x) , x — x}. The projec- 
tion of Tg onto the (x, u)-plane is the line x = x~. For the oscillator (|H|) 
this is shown , together with the projection of js, in Fig. |^. Therefore, 
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assuming u G [0, ir], 7 a intersects Tg at P and at P = (x, 7r — «, 0). In 
order to determine notice that 

(30) ^ = iMtan M = 
ax ri (x) 

r x (x) sin-u r. (x) 



n (x) y/l- S m 2 u n (x) yV (x) - a 2 

as it results directly from (|27|)i,2 and (|2|). In Qjjg) the choice of "+" 
(resp. "-") corresponds to u G [0, |[ (resp. w G]|,7r]). Hence, the 
possible position jump is described by the following formula: 

(31) 0-0 = A<f>(x,u) =±2|o| / r , _ da; 

A* ri (x) A/g 2 (x) - a 2 

where x* (see Fig. |5|) is the root of the equation 

(32) \ q { x )\ = \a\ 

and the sign + (resp. — ) corresponds to u G [vr/2, ir] (resp. u G 
[0, 7r/2]). In the case of the oscillator equation fl32|) becomes 

x = — — (1-x) 
\mc z roJ 

which has only one root x* in the interval [0, 1[. Relation (|3"0|) remains 
valid also for u G [— 7T, 0]. In this case the sign + (resp. -) corresponds 
to u G [-7T, — 7f/2] (resp. u G [-7r/2,0]). 

Since on a given characteristic only two points lie, P and P, belong- 
ing to the same energy level, a jump from P to P or vice-versa may 
happen only if one of these points is "in" while the other is "out". 
This occurs iff the function Xjj 1 (g), g (q,p) = being the equation of 
C {Si), takes opposite signs at points P, P, and we go to analyse when 
such is the case. 

It follows from ([El]) that for lagrangian ( |16|) 



(33) X? = v^ + v^ + *% A + g 2 A) 

oqi ag 2 r api ap 2 

Below we work with the local chart (r, <fi,y,a),y = (pi + p\) /m 2 c 2 , a = 
arctan — , on T* (M). In terms of these coordinates the Legendre map 
(pOD is given as follows: 

(34) y = i/j(r,x), a = 6 , 

with ip (r, x) = -^axL 2 x (r, x), and the expression (|33|) takes the form 

(35) 

„ rf , ,- ( d sin m d \ \pi\) _ / 9 sin w 9 \ 

X™ = Cv ^ cosm— + — + ^-L r 2cosm- — — 

\ or r o<p J vac \ ay ip oa J 
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FIGURE 6. Legendre map for the oscillator (* = r/rO) 



Due to 


©i 


hyper surface C (Si) is 


given 


by: 








g(r,y) = 


= 


i 


with 












(36) 




9 (r, y) 


4 

= 2 4 Xl 


(r)L 2 x 


(r) 



where x = x\ (r) is the function implicitly defined by equation E x (r, x) 
0, and L x (r) = L x (r, rci (r)). For the oscillator ( |j~4|) the Legendre map- 
ping £ is illustrated in Fig || In this case, as it is easy to see from ([34]) 
and (|36D , the image with respect to £ of a sufficiently small neighbour- 
hood of a point (q, v) G S\ is contained in the region {g > 0}. 
From ( j35|) we get: 

x h 1 (fiO = cosu( cy/xg r H L r g y J 

3^) and C (-S'i)'s equation follows: 



But from ([21 

4 

9r 



■m 2 c 4 



4L X 
m 2 c A 



X\E X + 2x\L xr 



m 2 c 4 



XiL x L x 



where, as before, / (r) means / (r, x\ (r)). Hence, we obtain: 

4 



(37) 



** 07) 



m 2 c 3 



y/xcosu \2xiL x L xr — L X L 1 



Since 2xL xr — L r = E r the restriction of ( j37|) to iS"i is: 

4 



(3 



^ (<?) 



m 2 c 3 



xi cos uL x E r 
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In the case of the oscillator we have L x < 0, E r > . Therefore (|38|) 

shows that X r jf l (g) is: 

1) positive in the region S± corresponding to u G [vr/2, 7r]U[— 7r, —tt/2]; 

2) negative in the region Sf corresponding to u G [— 7r/2, 7r/2]; 

3) null on the bidimensional surface W = Si H {u = ±7r/2}. 
Hence, keeping in mind what was said above about the range of C 

and applying the transition principle, we see that: 

1) If P G Si, then P is an in-point, while P is an out-point. There- 
fore, if the phase point, starting from outside S\, reaches it at P, its 
trajectory can be prolongated starting from P (jump from P to P). 

2)If PgSi, then P is an out-point, while P is an in-point. There- 
fore, if the phase point, starting from outside Si, reaches it at P, its 
trajectory can be prolongated starting from P (jump from P toP ). 

3)If P G W, i.e. if u = ±7r/2, then P = P. In this case the jump 
becomes infinitesimal and its direction is indicated by the hamilton- 
ian vector field. In fact £ (W) is described in terms of coordinates 
(r, 0, y, a) by equations 

9 (r, y) = 
cos u = 

But it follows from (|35f) that on W 



H (COSU) = pSin« = ± — , 

mcy/ip mc\/yj 

which is different from zero on Si. More precisely, it is positive for 
u = 7r/2 and negative for u = —vr/2; in both cases X^f 1 is directed 
toward the region Sj~ in which singular trajectories end. 
Finally note that, for the oscillator ([14]), 



X r «(pi) = ^L r ^0 on S 2 \S 



r 



Since C (S 2 ) coincides with the null section of T* (M), this shows that 
Xh is transversal to C {S-z). 

4.4. Phase trajectories of the oscillator. In this concluding sub- 
section we will study phase trajectories of the oscillator outside S. It 
will be shown that their behaviour depends strongly on their position 
with respect to the singular surface. Trajectories arriving at S have 
discontinuities, described by the transition principle, and together with 
others form a perfectly self-consistent dynamical model. 

In this subsection coordinates (r,x,u,<f)) are used. Let, as before 

„ . d ■ d , d , d 

Z l = r 7T + Vttt + x— + li- 
ar Off) ox ou 
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be the vector field on T (M) corresponding to the lagrangian L. Since 
/, E are first integrals, then 

Z L {I) = Z L {E) = , 

or, equivalently 





Therefore 

(39) r = -kI u E x , x = kI u E r , u = k {I r E x - I x E r ) 

On the other hand 

d ( nr< t\ ^i^i + ^2^2 i- 

T = ~dt I V g i +g 2 ) = = c^xcosu 

So that 

Cv^cosm Cy/xsinu c 2 

E X I U IE X 2rL x E x 

Notice also the relation: 

(41) I r E x — I x E r = — p= sinw (2xL x + rL r ) , 

C\/X 



which follows directly from (|i~8|) , (|21~D , (|19|) , and the relation 



(v^-^x) 



<9x Sy^x 

This way one gets the first three Euler-Lagrange equations: 
r = cy/xcosu 

rp 

(42) x = —c—^yfxcosu 

E x 

csinw 

u = -2^T x {2xLx + rLr) ' 

which form a closed subsystem of the whole system. 

Denote by Zl the projection of Zl onto the (x, r, u)-space. Then, 
solutions of ( f42|) are identified with trajectories of Zl. Due to obvious 
symmetry with respect to the (x, r)-plane it is sufficient to consider 
those of them for which u G [0, it] (i.e. counterclockwise motions around 
the center of the elastic force). 

The fourth Euler-Lagrange equation 



(43) = ^sinw 



can be found directly from : 

tan0 = 

9i 



tan <t> = — , 
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Figure 7. Energy level surfaces 



so that 

(44) 



6(t) = c I — r ^sinw(r) dr + const. 
r ( r ) 



with (x (t) ,r (t) ,u (£)) being a solution of (f42"|). 
Let 



(45) 
Obviously, 



2a, m = {-E = Amc 2 } n {I = fimcro} 



1 — 'A,/u 



A,/i 



xS 1 



where T\^ is the projection of H A M into (r, x, w)-space, while the circle 
S 1 corresponds to the cyclic coordinate (p. In their turn surfaces E\ ifl 
foliate the energy level 3-fold 

S A = {E = \mc 2 } 



In the case of oscillator ( |T4"D £a is not empty for A G [1, +00]. 

Let Sa be the projection of £a onto the (r, x, w)-space. Obviously 

s A = r A xS 1 , 

where Ta is the curve in (r, x)-plane given by equation E (r, x) = Xmc 2 
and S 1 is the circle corresponding to the cyclic coordinate u. Curves 
Ta are shown in Fig |7|. One can see that Ta intersects the projections 
Ci,C 2 , of 5*i, 5*2, respectively, as follows: i) at two different points 
Pi, P 2 , if A > 2; ii) at the single point Q, if A = 2; hi) nowhere, if 
1 < A < 2. Therefore £a intersects Si and S2'. 
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Figure 8. Phase portrait of the oscillator for a fixed 
value of the energy 



1. along two tori T t = £ A fi «% , i = 1, 2, if A > 2. These tori project 
onto Pj's and have (<j),u) as cyclic coordinates. 

2. along the circle 7 = Si fl^, with the cyclic coordinate 0, if A = 2. 

3. nowhere if A < 2. 

Therefore, T,\\S has three connected components, if A > 2 and is 
connected, if A < 2. 

In the case A > 2 the behaviour of phase trajectories depends strongly 
on the connected component of T,\\S they belong to. Due to (|44l) it 
is sufficient to study trajectories of (|42"D , i.e. connected components 
of curves T\^\S, with S being the projection of S onto the (r, x, u)- 
space. In Fig. |8| the projections into the (x, w)-plane of three different 
kinds of such trajectories contained in £a for a fixed A are shown. 
As before we limit ourselves to < u < n. The three vertical lines 
correspond to the projections Tj of tori Tj's (in fact they are circles, 
due to the ciclicity of coordinate u). Passing to further details, de- 
note by Xi = Xi (A) the constant value of x along T,. Connected 
components £^, E|, E| of E A \5 correspond to x < X\,X\ < x < x 2 
and X2 < x respectively. The trajectories belonging to H\ and E^ 
are discontinuous in the sense that they end on Si and then jump, 
according to the transition principle, at another point of Si. More 
exactly, in the situation shown in the figure, such a trajectory starts 
from a point of T± = T\ fl {7r/2 < u < 71} and reaches a point in 
T{ = Ti fl {0 < u < vr/2} in a finite time. Trajectories 71 and 72, 
whose projections 71 and 72 are shown in Fig. || correspond to the 
same value of A and of /1. Denote by P s tart = (%, r,n — u) G T x + and 
P end = (x, r, u) G Xf the common starting and ending points of 71 
and 72. When a phase point starts from (P s t ar t, 0o) = (x,r, ir — u, 4>q) 
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Figure 9. Jumping phase trajectories 




JUMP! 



Figure 10. Jumping oscillator 



and then goes along 71 (or, alternatively, 72) it arrives at the point 
(P en d, 0) = (x, r, u, (f) with given by ([H]) and then proceeds along 
the trajectory whose projection is 71 or, alternatively, 72, and so on. 

In Fig. [10] the situation in the configuration space M = {(qi,q 2 )} = 
{(r, 0)} is shown. The oscillating particle, starting from the point 
(r = r, <p = 0), moves along the projection 71 of 71 onto M until it 
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Figure 11. Regular precessions 

reaches the point A en( i = (r — f, (j) — 0] . From there it jumps to the 

point A star t= [r = r, = (f>), where <fi is given by (|31|). Then it splits 
into the two trajectories #i, 5 2 , both with initial velocity 



7T 



U 



contained in E|. These are 



Consider now the trajectories of 
regular closed trajectories winding around the center C = C (A) = 
(x = xc,u = vr/2), where xc = xc (A) is the zero of the equation 
2xL x (r en (x, A) , x) + rL r (r en (x, A) , x) = 0, and r = r en (x, A) is im- 
plicitly defined by the equation E(r,x) = Xmc 2 . The corresponding 
trajectories in the configuration space are shown in Fig. |11|. They are 
precessions around the force center. 

Note that both the discontinuous trajectories in T,\ and the regular 
ones in H 3 X near S2 tend to be parallel to it, so that this component of 
the singular hypersurface is never reached by the phase point, at least 
for trajectories with non-zero angular momentum. 
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